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In this lecture we present some interesting issues that arise when 
the dynamics of the charge carriers in the CuC>2 planes of the high 
temperature superconductors is considered. Based on the qualitative 
picture of doping, set by experiments and some previous calculations, 
we consider the strength of various inter and intra-cell charge transfer 
susceptibilities, the question of Coulomb screening and charge collective 
modes. The starting point is the usual p-d model extended by the long 
range Coulomb (LRC) interaction. Within this model it is possible to 
examine the case in which the LRC forces frustrate the electronic phase 
separation, the instability which is present in the model without an 
LRC interaction. While the static dielectric function in such systems 
is negative down to arbitrarily small wavevectors, the system is not 
unstable. We consider the dominant electronic charge susceptibilities 
and possible consequences for the lattice properties. 



1 Introduction 

A decade of experimental and theoretical research indeed showed that the physics 
of the high temperature copper-oxide superconductors is both complex and intrigu- 
ing. One major source of complexity lies in the fact that the electronic correlations 
in these materials are strong. The antiferromagnetism in undoped LaCu04 and 
YBa2Cu306 is clear evidence of that. A second source may be found in experiments 
that suggest that the electron-phonon interaction may be strong as well. As far as the 
effects of strong electron-electron interactions on the electronic properties of doped 



superconducting materials are concerned, two different viewpoints exist in general. 
From the first point of view, the spin correlations, although much weakened on dop- 
ing, still dominantly determine the carrier dynamics. From the second viewpoint 
the superconducting cuprates are doped charge transfer insulators, in which the ten- 
dency towards spin ordering may be an important, but still secondary, effect. This 
second approach focuses on the charge degrees of freedom and charge fluctuations, 
and the effects that they may have on superconductivity jj], ||, electronic properties 
in the metallic phase, and the crystal lattice dynamics Q. Some of these issues are 
the subject of this lecture. We consider dominant static charge susceptibilities, the 
collective modes and the dielectric properties for the three dimensional system of 
copper-oxide planes. The features that emerge are related to the properties of these 
materials observed in experiments. In particular, these are the anomalies found in 
the modes of oscillations of in-plane oxygen atoms, the apex oxygen position and 
movement and the incommensurate and commensurate deformations in the direction 
parallel to planes. In this lecture we consider the system that is close to the doped 
charge transfer insulator. However, some of the conclusions remain qualitatively 
valid away from this regime. This is in particular true for the foregoing discussion 
related to the oxygen-oxygen charge transfer modes and the discussion regarding 
the frustrated phase separated instability. 

2 The p-d model with long range Coulomb interaction 

The presence of CuC>2 planes is the mayor feature of all cuprate high temperature 
superconductors. The electronic hopping between the planes is relatively small. 
The 'insulating' layer between the planes may be composed by atoms of various 
elements. Another generic feature, as shown by many experiments (e.g., EELS, X 
ray absorption, photoemission), is that on doping the excess holes mostly populate 
the p Xj y — a orbitals of the in-plane oxygen atoms. These orbitals further hybridize 
with some, mainly d x 2_ y 2, orbitals on the copper atoms. Together with the fact 
that correlation effects are important for the in-plane dynamics of electrons, this 
constitutes the basis of the p-d model. 

The original p-d model [IJ consists of the tight-binding part and the part 
accounting for the short range Coulomb interaction, 



In the tight-binding part the orbital energy levels (e^, e p , A p d = e p — Ed) are spec- 
ified and various hybridization terms (copper-oxygen, to, oxygen-oxygen, t', etc.) 
included. Hsrc contains the terms describing the Coulomb repulsion on the cop- 
per, Ud-, and oxygen site, U p , as well as the terms describing the interaction between 
electrons on neighbouring sites (copper-oxygen V p d, the oxygen-oxygen V pp , etc.) 

In this lecture we will consider the extension of this model which includes the long 
range Coulomb forces. This extension is necessary in order to examine the capability 
of strongly correlated electrons to participate in the screening and to determine the 
effects that strong local forces have on electronic plasmon. The second reason for 
the inclusion of the long range Coulomb forces is to stabilise the system against the 
phase separation, the instability that occurs in the part of the parameter space of 
the original p-d model. 



H p d = H ot ba + Hsrc- 
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The introduction of long range Coulomb forces into the p-d model is relatively 
straightforward. However, a few points should be emphasized. First, as soon as 
long range Coulomb forces are considered, it is natural to extend the model from 
one plane to the three dimensional collection of parallel planes. The electron-electron 
interaction is described^ by 

H CO ui = - n a (R)V aa/ (R + a a - R' - a a >)n 3 (R') (2) 

R,R',a,a' 

where R = (Rn,R±) represents the cell index and a denotes the orbital of the 
atom positioned at a a with respect to the origin of the unit cell. The potential 
V a g(R + a a — ag), describes the interaction of electrons on different atoms in the 

unit cells separated by R. It is proportional to 1/ R + a a — ap at long distances 
and reflects the relative position and the charge distribution in the orbitals at short 
distances. The Fourier transform, V a $[k) = J2r V a p(R + a a — ap)exp(ikR) = 
V a g(k) exp(— ikag — ika a ) behaves like 1/k 2 at long wavelengths, while the short 
wave dependence reflects the structure and strength of the local forces. This may 
be emphasized by writing V a g(k) in the formP] 

-» 47re 2 -» 
VaA k ) = -71 7fZ + CofiW (3) 

where a 2 d± is the volume of the unit cell (a is the lattice constant of the Cu02 
plane and d± is the distance between neighboring planes). The high frequency 
atomic screening is included through e^. The terms C a p(k) have finite values^] at 
k = 0. 

As usual, the k = term of the 1/k 2 part of the Hamiltonian cancels out because 
of the neutrality of the crystal. Further on, the 'passive' orbitals on the ions with 
fixed valence may be excluded from the electronic Hamiltonian, leaving only the 
active orbitals in the copper oxide planes in the model^j 

The Coulomb part of the p-d model, with the long range forces included, finally 
reads as 



n coul / j 2 

k£B.Z.,a,f3 



(1 -%o)^T + ^) 



n a (k)n B (-k) (4) 



As in the original p-d model, we neglect the Fock terms here. Physically, this is correct for the 
long range, but not necessary correct for the short range Coulomb forces. 

2 A somewhat more complicated, but less suggestive form should be preferable for large fc, 
close to the zone boundary. For example, an additional exp(— k 2 a 2 ) factor may be introduced 
in order to make the cutoff at the zone boundary soft. Also, it is possible to replace k 2 by 
Si=a; y z — cos(fc;d;)] /df in order to have zero derivatives at the zone boundaries. 

3 It is an instructive exercise to reformulate the calculation of the Madelung energy of some pure 
ionic crystal after the outlined procedure. The contribution of the 1/k 2 term then explicitly cancels 
out due to the neutrality of the crystal, ^ n a = 0, and the Madelung energy per unit cell acquires 
the form E M /N = (1/2) C a/3 (O)n a n0. 

4 The terms neglected should be kept in mind as the source of a shift, presumably small, of 
the energy levels of active orbitals with doping. More important, these terms reappear as the 
electron-phonon coupling terms when lattice movement is considered. 
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where n a (k) oc J2r n a(R) ex P( — ikR — ika a ) and the dimensionless parameter w = 
e 2 /<i_i_to£oo> measuring the strength of the long range Coulomb interaction is intro- 
duced. The approximate value for w expected for cuprate superconductors may be 
estimated by taking d± ~ 12A, Soo ~ 3, to rs 1.3eV. This gives w ~ 0.25. 

In the text that follows we will consider the situations with and without the 
1/k 2 term in the Hamiltonian. In the latter case, C a p(k) corresponds to Ud, V p d, 
Vp P usually used in the literature. Physically, they represent the effective, screened 
electron-electron interaction. In the case where 1/k 2 is explicitly included we con- 
tinue to use the notation Ud, V p d, V pp , etc., although the physical meaning of C a/ g's 
in that case is somewhat different. They represent the local field corrections to the 
1/k 2 interaction, presumably large for short distances. 



3 Qualitative features of the electronic spectrum of doped 
charge transfer insulator 

Various estimates of the local Coulomb terms in the p-d model emphasize the large 
repulsion on the copper site, Ud ~ 10 eV. The limit Ud — > oo is frequently considered 
in the literature, and we will work in that limit as well. The rid < 1 constraint forced 
by infinite Ud is often dealt with by introducing the auxiliary (slave) boson field b 
and the auxiliary field fermion / in terms of which the real electron operator on the 
copper site becomes the composite object, dg = bj-ffi. The next step is usually to 
consider the mean field approximation for the boson field (the saddle point approx- 
imation in the partition function/path integral language). Usually, the corrections 
beyond the mean field (the Gaussian fluctuations around the saddle point) are also 
considered. While rigorous justification for this procedure is usually sought in the 
1/N expansion of the generalised model (i.e. N spin components instead of two), the 
physical appeal comes from the fact that it grasps some major physical features. In 
particular, this applies for related models Q|, that may be more accurately investi- 
gated by other means (for instance, the appearance of the Abrikosov-Suhl resonance 
in the Anderson model). For the p-d model this procedure qualitatively reproduces 
some basic experimental facts in cuprates as are the appearance of the in-gap res- 
onance H and the large Fermi surface |?J on doping. A pedagogical survey of the 
method and some of its results may be found in ref. ||. 

What basically happens at the mean field level H for the slave boson in the p-d 
model with large Ud is: a) the boson field b acquires a static component Bq (with 
rid = 1 — \Bq\ 2 < 1) and an oscillatory component at the frequency A; b) A also 
represents the shift of the copper orbital energy on changing from original electron 
to auxiliary fermion fields, £/ = £4 + A; the copper-oxygen hybridization in the 
auxiliary fermion Hamiltonian is reduced with respect to the original hybridization, 

The band that occurs at e ~ Et (/-band) in the auxiliary fermion density of 
states becomes the in-gap resonance when the real electron spectrum is considered. 
Quite simply, the majority of the spectral weight for the copper site that lies in the 



auxiliary fermion /-band is transfered back [10, |4| to original energy Ed = et — A on 
calculating the auxiliary fermion - auxiliary boson convolution in the real electron 
Green function (-«) (T(b^ f a )(t)(pb)(Q)\. The spectra are illustrated in Fig. [I]. 



The Fermi energy stays in the resonance. The Fermi surface states are predom- 
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Figure 1: a) The density of states in doped p-d system, as obtained at the mean field level, 
b) The density of states for the auxiliary fermions. The labels (p)[t](d) indicate the total 
weight of the particular part of the spectrum [t] and the shares (p) and (d) of the oxygen 
and the copper orbitals, respectively. The resonance at e rj £/ in (a) and the corresponding 
band in (b) are approximately half filed for small doping level 5. of the oxygen (p) and the 
copper (d). 

inantly made of oxygen orbitals, especially in the doped charge transfer insulator 
(CTI) regime (A p d ^> to). In that regime the mean field value Bq vanishes when the 
doping 5 (measuring the hole concentration with respect to the level of one hole per 
cell) goes to zero, Bq oc 5. The renormalized auxiliary fermion bandwidth W oc Bq, 
being also the effective width of the resonance, and the number of states in the 
resonance are also proportional to doping level 5. 

It will be useful to keep this qualitative picture^] in mind when the results for 
the charge susceptibilities are discussed. 



4 Charge transfer susceptibilities and Coulomb 
screening 

The results for the charge susceptibilities that we present are calculated on taking 
into account the first fluctuation correction beyond the mean field approximation. 
While for large Ud the slave boson approach is used, for other, smaller Coulomb 
terms we use the Hartree^j approximation |8|, 12] as the mean field approximation. 
The fluctuation correction for the long range Coulomb interaction corresponds to 
the usual random phase approximation (RPA). Fig. ^ reflects the simple structure 
of our calculations. The charge transfers that we will consider may be expressed 
in terms of the copper and the oxygen site charges rid, n x , n y . The combinations 
reflecting different charge transfer symmetries and different physics are n p = n x + n y 



(the charge on oxygen sites), 



n 



pp 



n. 



- n, 



(the charge transfer between inequivalent 



oxygen sites), n pd 



rid (the copper-oxygen charge transfer) and n n = n p + rid 



(the total charge in the cell). The susceptibilities Xa = {(n a n a )), ol = d,p,pp,pd,n 

5 This picture may be refined by calculating the auxiliary field propagators beyond the saddle 
point level ^, [Hj and by calculating the electronic spectrum beyond the decoupling approxima- 
tion, leading to simple convolution. 

6 The Hartree approximation and the RPA fluctuation correction does not seem to provide an 
adequate treatment for the short range Coulomb interaction terms as is the V p d term. However, we 
will really consider only some qualitative features that these interaction bring in. These features 
are expected to remain present even if some more adequate treatment like the full Hartree-Fock 
approximation is used. 
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Figure 2: The corrections beyond the saddle point and Hartree approximations. 

will be considered. We will distinguish the susceptibilities x° of the tight binding 
fermions with the mean field renormalized band parameters; the susceptibilities x S 
calculated for the model with short range forces; the susceptibilities x calculated for 
the model in which the long range Coulomb forces are included as well. 

In Figs. H various static charge transfer susceptibilities are shown as a function 
of the wave vector k = (k, 0, 0) (hereafter we will use the bare copper-oxygen hy- 
bridization energy t as the energy unit and the Cu02 plane lattice constant a as 
the unit distance in order to simplify the notation) . From these figures the following 
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Figure 3: Static charge susceptibilities for A p d — 4,6 — 0.1 as functions of the wavevector, 
k = (k, 0,0): Shown susceptibilities reflect the transfers of a) the total cell charge b) total 
oxygen charge in the cell c) the copper charge d) the intracell oxygen-oxygen charge transfer. 

may be noted: 

a) The inclusion of the fluctuation of the slave boson field, corresponding to Ud = 
oo, suppresses the charge fluctuation on the copper sites for all wavevectors: 
Xd ^ Xd ■ The copper-oxygen charge transfer susceptibility is also diminished. 
It can be also seen that the Fermi surface enhancement]] seen in \d a ^ k x ~ 0.1-7T 

7 Several types of the Fermi surface enhancements of the static susceptibilities in the two dimen- 



G 



is lost in Xd- 

b) The fluctuation of the oxygen charge n p is not so much affected by the auxiliary 
boson fluctuations. However, the Fermi surface effect disappears in Xp as well. 

c) The inclusion of the long range Coulomb forces further suppresses the total 
charge fluctuations at small wavevector, Xn "C Xn- While this is expected to 
happen, it is important to realize that the charge fluctuation on the oxygen 
sites are mostly affected. This implies that the holes in the oxygen orbitals 
dominantly participate in the Coulomb screening and that their dynamics 
determines the dielectric properties of the system. This issue will become 
more clear when we consider the charge fluctuation spectra. 

d) The pp susceptibility related to the charge transfer of the between inequivalent 
oxygen atoms (x and y) in the unit cell is unaffected both by the long range 
Coulomb interaction and by the slave boson fluctuations. The reason lies in 
the particular symmetry of the pp charge transfer. In particular, the Fermi 
surface effects are pronounced in Xp P an d Xppi m contrast to the their absence 
in the susceptibilities Xp an d Xp related to the total oxygen charge. In fact, all 
van Hove and nesting enhancements present in the tight binding's x° remain 
in Xpp whereas they are suppressed in all other channels by the large Ud- This, 
in particular, means that the k = charge transfer transitions |l3|, and 
k = (-7T, 7r) charge density wave should be searched for mainly in the pp channel. 
The doping level particularly favorable is the one for which the Fermi surface 
touches the zone boundary. 

The dependence of the static dielectric function on the wavevector is shown in 
Fig. ||. The behaviour l/e(k) cx k 2 — > as k — > corresponds to the full metallic 
screening. More important, it should be noted that the screening distance, expressed 
as l/k s in Thomas-Fermi language, where l/s(k) = k 2 /{k 2 + k 2 ), is of the order of 
the lattice spacing, in spite of the fact that the concentration of doped holes is 
small (average distance between them is approximately three lattice constants for 
5 = 0.1). The doping dependence of the screening length shows no drastic changes 
with doping. This reflects the fact that density of states of the in-gap resonance (i.e., 
(number-of-states) /(resonance- width)) does not change very much with doping, even 
in the doped CTI regime. 

5 Dynamic correlation functions and collective modes 

More detailed information on the charge dynamics are obtained on considering the 
dynamic charge correlation functions. They further clarify which orbitals participate 
in various intra and inter-cell charge transfers and how they are affected by short and 
long range Coulomb forces. Some aspects of the the spectrum of charge fluctuations 
for the system with large Ud was considered by Kotliar, Castellani and coworkers 
H], [ilf . They discussed the appearance of the high frequency exciton collective mode 
{oj = ujpd) as well as the zero sound mode in the large Ud p-d model. The spectrum 

sional tight binding model occur at k « and k « (±7r, ±7r): van Hove effects, as well as the nesting 
or the sliding Fermi surface effects, if some flat parts of the Fermi surface exist. The degeneracy of 
these effects for the square Fermi surface is lifted for some more complicated Fermi surface shapes. 
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Figure 4: The figure on the left shows the dielectric function for (a) V p d — and (b) 
V p d = 1.25. The second case (as discussed in section 6.) is rather close to the border of the 
'normal' (V p d < 1.15) and FEPSI (V p d > 1.15) situation. The figure on the right shows the 
total charge susceptibility in the unstable V p d = 1.25 system without the long range forces, 
X^(0) < 0. Other parameters are: A p d = 4, 6 = 0.1, V pp = 0. 



of the oxygen charge fluctuations shown in Fig. || the appearance of these modes 
in the spectrum together with the intraband and interband continua. Fig. |5|) 
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Figure 5: Left: The spectrum of the fluctuation of the charge in the oxygen orbitals, 
n p = n x + n y . The full line correspond to the case without long range interaction. The inter- 
cell and the intra-cell (Cu-O) charge fluctuations, corresponding to the zero sound mode and 
the Cu-0 exciton mode, and may be distinguished. The long range Coulomb interaction 
changes the zero sound to the intraband plasmon mode, pushing its frequency to finite 
frequency at low wavelengths. Right: In the intra-cell oxygen-oxygen charge fluctuation 
spectrum, n pp — n x — n y , there are no collective modes present. The strength of the 
interband contribution in \pp is big relative to those in other spectra. Parameters as before. 
V pd = 0, k/ir= (0.11,0,0). 

also shows the effect of the long range Coulomb interaction on the spectra. It is 
interesting to note that the frequency of the exciton mode is not at all affected 
at long wavelengths. This shows the quadrupolar character of the corresponding 
charge excitations. On the other hand, zero sound mode changes to the intraband 
plasmon mode with the dispersion uj\{k) = a;f p fc|/(fc| + k^) usual for quasi-2d 
systems with negligible inter-plane hopping. The frequency (as well as the strength 
of the corresponding pole in the density-density correlation function) vanishes for 
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Figure 6: Left: The spectrum of charge fluctuations on the copper site. The Cu-0 exciton 
mode and the interband excitations are pronounced. Right: The total charge fluctuation 
spectrum as seen through l/s(k,u>). The intraband plasmon and the intra-cell dipolar 
excitation, which lies within the interband continuum, dominate the spectrum. 



the wave vector directed perpendicular to the planes. The frequency uj\ v of the long 
wavelength in-plane oscillations increases on increasing the strength w of the long 
range Coulomb interaction, but always stays in the gap. The dependence on w may 
be understood through a simple equation for the Coulomb collective modes in the 
two (infinitely narrow) band systems, 

uj z ui - 

Here A p f = Ef — e p measures the renormalised gap while Vl\ and S7| are products of 
the Coulomb factor w oc e 2 and the intraband and the interband oscillator strengths, 
respectively. reflecting the intraband charge transfer oscillator strength, assumes 
particularly simple doping dependence Sl\ oc 5 in the doped CTI regime. This is 
indeed expected from the picture that we have developed until now. 

The second solution of the equation (corresponding to the third collective mode in 
the calculation) represents the intracell dipolar mode. It starts from finite frequency 
lu ~ A p f for small w and becomes the 'big plasmon' as w is further increased. 
However, in our calculations the substantial w > 1 is required to push the frequency 
of this mode above the interband excitation continuum. 

In more detailed inspection of the spectra ^{k^uj) = lmx a (k,uj) of various 
charge transfer modes, we find that the intraband part of the spectrum is suppressed 
in all channels, except from the pp modej^j The zero sound/plasmon mode dominates 
in the low frequency part of spectra of the total charge fluctuation x^(/c,a;) and the 
oxygen charge fluctuation Xp(k, w). The strength of this mode in the channels related 
to long wavelength charge fluctuations on the copper site is rather small. These 
spectra are dominated by the interband excitations and the quadrupolar exciton 
mode. The later is absent in the total charge fluctuations at long wavelengths. The 
imaginary part of the inverse dielectric function^] is shown in Fig. ||. 

8 None of the collective modes appears in X^p(k,tj) for k = 0. For finite k their strength is 
negligible. 

9 We extract the 'macroscopic' dielectric function from the matrix of charge susceptibilities by 
considering the coupling of the system to the external potential. 
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6 Frustrated electronic phase separation instability (FEPSI) 



The dynamics and the dielectric properties of the system with frustrated electronic 
phase separation instability is even more interesting to consider. Experimentally, 
the proximity of the phase separation instability (PSI) and superconductivity in 
cuprates and related materials is rather well documented The possibility for 

PSI in the model with strong electron-electron interaction was addressed by several 
authors. In particular, V.J. Emery and S.A. Kivelson considered the case of low 
doped antiferromagnet and discussed the possible consequences of FEPSI for the 



normal state properties and superconductivity in cuprates |17|| . 

The phase separation instability within the p-d model with short range forces 
occurs as the copper-oxygen Coulomb term V p d is increased in the hole doped system, 
as pointed out in ref. |l2| . However, the enhancement of the copper-oxygen charge 
transfer by Vpd, often emphasized by some authors fil 0|j is n °t an important issue 
here.00 

The main reason for a system with strong correlations to choose the phase sepa- 
rated phase is to diminish its kinetic energy. The kinetic energy diminishes on doping 
since the frustration for the electron hopping caused by strong coulomb interaction 
becomes less effective]^]. On approaching the point of thermodynamic instability 
the derivative (dfi/dn) = (d(j,/d5) (the rate of change of the chemical potential with 
doping) approaches zero and becomes negative in the unstable phase. Also, the 
susceptibility Xn(k = 0,cj = 0) (as well as all other Xa s ) diverges at this point, 
1/Xn{k = 0,u; = 0) = (dfi/dn). An example of l/x„(/s) in the unstable system is 
shown in Fig. 

Of course, once the long range Coulomb interaction is introduced into the model, 
the true phase separated state is not likely to appear because of the huge costs in the 
Coulomb energy. Speculating on this issue, various authors stop at this point and 
suggest that the charge density wave (CDW) state will replace the phase separated 
state (i.e. that the instability will show as the divergence of Xn at some finite 
wavevector instead of at k = 0). However, the pure electronic CDW (assuming that 
the lattice is too rigid to participate in the formation of the CDW) corresponds to 
unrealistically weak Coulomb forces for the cuprate superconductors, id « 1. We 
find that the stable FEPSI state with l/x a (k) > and homogeneous electronic 
density is more probable. However, the static dielectric function, 

i i 

(6) 



e(k) l + (47™;/fc 2 )x£(fc||) 

is negative in FEPSI systems at long wavelengths since Xn(^)' which accounts for 
the local forces, turns negative, l/x^(0) oc (dfi/dn) < 0, as shown in Fig. |j. 

10 The value of V p d for which the enhancement of the copper-oxygen charge fluctuation takes place 
(possibly related to the softening of the corresponding exciton mode) is much bigger than the value 
required for PSI. For example, in our calculations for A p d = 4to PSI occurs already at V p d = 1.15to, 
while a value approximately two times bigger is required for the cooper-oxygen charge transfer 
instability to happen. 

11 The p-d model in the parameter range (Ud — A p( j) <C A p< j <C Ud exhibits the phase separation 
instability Jl(J even for Vpd = 0. 

12 V.J. Emery originally discussed and repeatedly emphasized this issues in the framework low 
doped t-J model. 
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The figure suggests the following formula to describe the dependence of the 
dielectric function on the wavevector, 



e(k) ~ 1-6 \k 2 + k 2 sl 



k 2 b 
k 2 + k 2 



(7) 



s2 , 



with l/k S 2 <C k S 2 {\/k S 2 is at least of the order of few lattice constants in the CuC>2 
plane) and b < 1. The formula resembles the Thomas- Fermi formula, but has two 
characteristic wavevectors instead of one. In order to get a qualitative feeling for the 
system, one may try to play with a three dimensional toy model with the dielectric 
function of the form (|7|) . The feature that readily emerges is the overscreening of the 
test charge at distances l/k s \ followed by the complete screening at distances beyond 
l/k S 2- This implies the attraction of two equal test charges down to distances l/k S 2- 
Also, the limit 1/e — > as k — > accounts for the total metallic screening. It is easy 
to see that a piece of the material exposed to the static electric field behaves just 
like an ordinary metal, developing a finite surface charge in order to ensure E = 
in the bulk. 

Returning to our quasi-two dimensional FEPSI metal we examine its static 
charge susceptibility. It is given in a form that may be easily analysed (see again 
Fig. 1), 

1 4ttw 1 . . 

+ — ■ (8) 



Xn(k) k 2 ± + kj Xn( k \\) 

The most pronounced features (see Fig. |7|) are the maximum at finite fen = k c (not 
related to some particular Fermi surface wave vector - k c decreases with decreasing 
w) for k± = 0, and the increase of Xn towards the k± = ir/d_i_ zone boundary. It 
may be noted here that the divergence of Xn{k\\ = k c ) for sufficiently small w is the 
usually mentioned CDW instability which comes as the alternative to the FEPSI 
homogeneous in the Cu02 plane. This instability corresponds to the softening of 
the intraband plasmon branch at fey = k c . However, the Cu-0 exciton mode is not 
affected at this point. 
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Figure 7: Left: The dependence of the total charge susceptibility on k±, k±/ir — 
0(a), 0.1, .., 0.6(g). Right: The possible effect of large inter-plane charge transfer suscep- 
tibility on apex oxygen atoms in the systems with CuC>2 bilayers, like YBa2Cu307_a;. 
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7 Lattice properties and the stability of FEPSI systems 



Finally, we may turn to the question of the stability of the FEPSI system, consid- 
ering in particular the influence of the electronic gas on the ionic lattice (which was 
considered infinitely rigid up to now). First, the question of total compressibility 
of the FEPSI system may be addressed. It may be shown in a number of ways 
that the short range forces between ions provide the stability of the system. One 
way is to consider the longitudinal sound velocity and the corresponding frequency 
uj 2 {k) = c 2 k 2 = Cpk 2 + £lp tian / 's e i{k) consisting of the contribution of the short range 
forces between ions and the long range forces screened by electrons. For FEPSI 
systems the second term is negative. In terms of the total bulk modulus B which 
determines the sound velocity c 2 = B/p, the contribution of the second term equals 
n 2 (dfi/dn) < 0. This, however, does not jeopardize the stability of the system, since 
size of the the electronic contribution is two orders of magnitude smaller than the 
total bulk modulus measured for cupratesj^] B ~ 10~ 12 dyn/cm 2 . The reason for 
small absolute value of the electronic contribution may be ultimately traced back to 
the small density of the electronic gasj^ 

More interesting than the compression of the whole system are the phonon modes 
that couple directly to the charge fluctuations that our model calculations distin- 
guishes as relatively big. In summary^, these are: a) the in-cell oxygen-oxygen 
charge fluctuations b) the k± = tt inter-plane charge transfer (this becomes k± = 
when systems with Cu02 bilayers are considered) c) the incommensurate charge 
density fluctuations inside the Cu02 planes. These modes may be expected to be- 
come soft Cj 2 h = oj 2 h — g 2 pXei when the electron-phonon coupling g ep is substantial 
or lattice instabilities and/or anharmonicities may result. 

In that respect it is interesting to note that, indeed, in several copper oxide 
superconductors anomalies in the lattice modes coupled to the intracell oxygen- 
oxygen charge transfer are observed. For example, in T^P^CaC^Og in the normal 
state the dominant mode for the in-plane oxygen motion is one that couples to 
n p p. This kind of oscillation becomes suppressed [19] in the superconducting phase. 
Similar competition of the superconductivity and lattice deformation that couples 
to n p p is found (2^, ||] in the LTT phase of La2_ x Ba z Cu04. 

The discussion on the experimental evidence for the in-plane incommensurate 
CDW in related, but non-superconducting compounds, may be found in ref. 17]. 

Here we would like to point out one possible sign of the FEPSI situation in 

13 The correct order of magnitude for the crystal compressibility may by be obtained from simple 
Madelung calculations. 

14 In the jellium language, r s ~ 8, if we count one hole per unit cell; counting only doped holes 
we get r s ~ 16 for 5 — 0.1. At this point, it seems interesting to note that the theoretical search for 
the high temperature superconductivity before the discovery of cuprate superconductors pointed, 
indeed, towards the low density metals M]. One of the basic properties of these systems is the 
appearance of the negative static dielectric function. That negative static dielectric function is 
the feature desirable for the electronic mechanism of superconductivity was rather pedagogically 
explained by Littlewood in ref. ^j. However, as emphasized by Kirzhnits at al. Jlq| , this is not 
the sufficient condition for the superconductivity to occur - the interaction between test charges is 
not the one that enters the gap equation - the local field corrections and correlation effects that 
themselves lead to negative static dielectric function should be taken care of, as well. 

15 The possibility of soft Cu-O charge transfer exciton is out of focus of this lecture. However, the 
formalism set here is perfectly suitable for its consideration since, as already explained, one almost 
inevitably works 'deeply' inside the FEPSI phase. 
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YBa2Cu307-x) the compound which contains two Cu02 layers per unit cell. From 
our discussion we expect to find anomalies for ion movements that couple to the 
charge transfer between these layers. The IR- active mode of the apex oxygen atoms 
is a natural candidate for such coupling. The coupling, being strong enough, results 
in the situation with two equivalent, minimum energy positions for the apex oxygen 
atoms (see Fig. [?| for illustration). Experimentally, this picture, supported later on 
by some other measurements, first emerged from the analysis of the EXAFS [21] in 
YBasCusCV, 



8 Conclusion 

The strong short range Coulomb forces substantially complicate the charge dynamics 
in doped copper-oxide planes. It seems difficult to envisage some effective band 
picture which would simultaneously account for the static (k s ) and dynamic (u p ) 
screening as well as some Fermi surface effects that we pointed out. Also, the strong 
Coulomb forces may tend to cause electronic phase separation, the vicinity of which 
was pointed out by several authors as the possible source of superconductive pairing 
and the anomalous normal state properties. At this point the introduction of the long 
range Coulomb forces into the model seems crucial, what we did for the p-d model. 
We showed that, while the long range Coulomb forces suppresses the instability, 
the negative static dielectric function and, more important, rather big total charge 
susceptibilities in some parts of the A;-space remain as the characteristic signs of 
FEPSI systems. Some of these signs were found in cuprate superconductors and 
related materials. While the lattice instabilities which may occur in order to exploit 
the large electronic susceptibilities probably do not contribute to the appearance of 
superconductivity, the effective electronic interaction which drives FEPSI may be 
favorable in that sense. 
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